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Abstract 

We analyse the CP asymmetry for Bg — > 77 in the two Higgs doublet model with 
tree level flavor changing currents (model III) and three Higgs doublet model with O2 
symmetry in the Higgs sector, including O7 type long distance effects. Further, we study 
the dependencies of the branching ratio Br{Bs — > 77) and the ratio of CP-even and 
CP-odd amplitude squares, R = \A'^\^ /\A~\'^ , on the CP parameter sinO. We found 
that, there is a weak CP asymmetry, at the order of 10"^. Besides, the branching ratio 

Br(Bs — > 77), and also R ratio, is not sensitive to the CP parameter for | - o'" | < 1. 



*E-mail address: mugeboz@thepl.hun.edu.tr 
tE-mail address: eiltan@heraklit.physics.metu.edu.tr 



1 Introduction 



One of the most important classes of decays are rare B-meson decays which are induced by flavor 
changing neutral currents (FCNC) at loop level in the Standard Model (SM). Therefore, one can 
obtain quantitative information for the SM parameters , such as Cabbibo-Kobayashi-Maskawa 
(CKM) matrix elements, leptonic decay constants, CP ratio, etc. and the measurement of 
physical quantities like branching ratio (Br), CP asymmetry (Acp),--- etc, gives important 
clues about the model under consideration. These decays are also sensitive to the physics beyond 
the SM, such as two Higgs Doublet model (2HDM), Minimal Supersymmetric extension of the 
SM (MSSM) 0. 

Bs ^ 77 decay, which is an example of the rare B meson decays, is an important candidate 
to test the theoretical models and to construct new models in the framework of the planned 
experiments at the upcoming KEK and SLAC-B factories and existing hadronic accelerators. 
Bs — ^ 77 decay, induced by the process b — » 577 in the quark level, has been studied in the 
SM 0-0 and 2HDM without QCD corrections. Since the QCD corrections to the inclusive 
decay 6 — > S7 are large (see - [|1^ and references therein), it is expected that they are also 
large in the inclusive b s'j'j decay and, therefore, in the exclusive B^ — > 77 decay. With the 
addition of the leading logarithmic (LLog) QCD corrections, the analysis has been repeated 



in the SM [[TTl-lTg, 2HDM [0], MSSM |T5| , and the strong sensitivity to these corrections 
was obtained. Recently, Bs — > 77 decay has been calculated in the framework of model III 
||T6| with the addition of the LLog QCD corrections and the current upper limit of the Br 
ratio Br{Bs —>■ 77) < 1.48 .10^"^ |]T^ was theoretically obtained for larger values of the Yukawa 
coupling ^N,bb and the ratio \rtb\ = ||^| >> 1 (see Appendix C for definitions). 

In the present work, we study Bg 77 decay in 2HDM (model III) and three Higgs 
doublet model with O2 symmetry in the Higgs sector {3HDM{02)) fl^, with complex Yukawa 
couplings. In our calculations, we take into account the perturbative QCD corrections in the 
LLog approximation by following a method based on heavy quark effective theory (HQET) 



for the bound state of the Bs meson [11 and we also consider the long-distance effects due 



to the process Bs ^ 4>'y ^ 77, [11 1. Since the complex Yukawa couplings are chosen in the 



calculations, in addition to the CKM matrix elements, there is a new source for CP violation. 
Using this new source, we obtain Acp in model HI at the order of 10~^, which is a small 
effect. Furthermore, we calculate this quantity in the 3HDM{02) and find that there is a 
small decrease as compared to the former one for sin 6 < 0.5. Finally, the calculations of Br 
and the CP ratio i? = [^[ show that these physical quantities are not sensitive to the models 
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under consideration for \rtb\ = \ < 1 (see Section 3). 

The paper is organized as follows: In section 2, we present the LLog QCD corrected am- 
plitude for the exclusive decay Bg — > 77 and the CP-even and CP-odd amplitudes 
in a HQET inspired approach. Then, we calculate Acp, assuming that the Yukawa couplings 
are complex in general and also derive Br. Finally, these calculations are repeated in the 
3HDM{02) model with the redefinition of the Yukawa combination Xg (see eq. (^)). Section 
3 is devoted to our analysis of the physical quantities under consideration and our conclusions. 
In the Appendix, we give a brief explanation about the model III and 3HDM{02) which we 
study. Further, we present the operator basis and the Wilson coefficients responsible for the 
inclusive b —>■ 577 decay in the model III. Finally, we give the explicit forms of some functions 
appearing in the Wilson coefficients. 

2 Leading logarithmic improved short-distance contri- 
butions in the model III for the decay Bg 77 

The LLog corrected effective Hamiltonian in the model III (see Appendix A) for the exclusive 
Bs 77 decay is 

where the Oj, O'^ are operators given in eqs. (^1]), (^21), and Cj, C- are the Wilson coefficients 
renormalized at the scale fi. The Wilson coefficients Cj and C- can be calculated perturbatively 
and their explicit forms of at mw level are presented in Appendix B. The effective Hamiltonian 
(0) is obtained by integrating out the heavy degrees of freedom, i.e , t quark, W^, H^, Hi, and 
H2 bosons in the present case where denote charged. Hi and H2 denote neutral Higgs 
bosons. Further, QCD corrections are done through matching the full theory with the effective 
low energy theory at the high scale /i = mw and evaluating the Wilson coefficients from mw 
down to the lower scale fi ~ 0{mf,). In our case, we choose the higher scale as /i = mw 
since the evaluation from the scale /i = mH± to /i = mw gives negligible contribution to the 
Wilson coefficients (~ 5%) since the charged Higgs boson is heavy due to the current theoretical 
restrictions (see fl^, ) 



The decay amplitude for the Bg 77 decay is obtained by sandwiching the effective 
Hamiltonian between the B^ and two photon states, i.e. < Ps|7ie//|77 >, and the matrix 
element can be written in terms of two Lorentz structures 0-0, [pl]] -|[13||: 



A{Bs 77) = J-^,^'^^ + lA^J^^,^^" , (2) 



where T^^y = ^e^^a/s^"'^ and {A~) is the CP-even (CP-odd) amphtude. Denoting CP-even 
(CP-odd) amphtude coming from the first operator set eq. (^Tj) as Af (Ai) and the primed 
operator set eq. ( ^2]) as A2 (^2 )> a HQET inspired approach, we get, 



At 
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Q eff I eff 
y rUf^ + TTLs 



(-mb J (m^) nisJ {ms))D{^) - mj {mc)E{fi) ' 



V2: 



fB^t 77 



3mB,As{mB, - 



(?_C7^^^(/x)-5:Qj/(m,)C,(^) 



and. 



+ 

At 
A, 



(3) 



TT 



1 



c;^^^(/i) 



(4) 



3mB,As{mB, - A 

However, we do not take the CP-even and CP-odd amphtudes corresponding to the primed 
operator set since their contribution is smaU as compared to former ones (see [^, |20|). In eqs. 
(I) and (0), Qf, = I for g = u, c and Qg = for q = d, s, b. Here, we have used the unitarity 
of the CKM-matrix J2i=u ct^H^ib = 0, and also the contribution due to V*gVub <^ VZ^tb = 



is neglected. The function g_ is defined as ITT 

9- = ^B, 



b "T Ill's ) ) ■ (5) 

The parameter A^ enters in eqs. (^ and (|) through the bound state kinematics (for details 
np. In the expression (^, the LLog QCD corrected Wilson coefficients Ci...io(/^) mi - 



see 



13| are, 

C,(/i) = {CM - C,{ii))N, + C4(/i) - C6(/i) , 

{CM + CM - cM)Nc + C2(/i) + CM - CM 
CM = {CM) + cM)){Nc + 1) - N,CM - CM 
cM + cMNc, 
cM + cMNc 

and the effective coefficient C^" {^), defined in the NDR scheme, is 

C^,"{^,) = Cr^'^{^^) + Q, (C|™^(/i) + C6^™^(/.)) , 

+ g„(^C?o™^(^) + Ar.!^C,^™^(;.)) 
nib rrib 



CM) 

CM 
Cs{^^) 
D{fi) 
E(/i) 



(6) 



(7) 



where A^^^ is the number of colours {Nc = 3 for QCD). The functions I {nig), J{mq) and A{mg) 
come from the irreducible diagrams with an internal q type quark propagating and their explicit 
forms are given in Appendix D. In our numerical analysis we used the input values given in 
Table (|lD. 



I—' Q T" Q m 

1 (XL ailic Lei 


Value 


rric 


LA V j 


mb 


4.0 (^^oreV j 




1 OQ 


At 


U.U4 




ZL DQ 1 n-13 ('Pf^A/'i 
^t.Uy ■ iU V J 


IBs 


U.Z, I VJC V I 




5.369 (GeV) 


mt 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


A (5) 


0.214 (GeV) 


Ois{mz) 


0.117 


As 


0.12 (GeV^) 


Ai 


-0.29 (GeV^) 


As 


590 (MeV) 


A 


500 (MeV) 



Table 1: Values of the input parameters used in the numerical calculations unless otherwise 
specified. 



At this stage, we will calculate the CP violating asymmetry for the given process. The 
possible sources of such effects are the complex Yukawa couplings in the model III. Here, we 
neglect all the Yukawa couplings except ^^^^^ and i^,bb (see Section 3). Using the expressions 
for the decay amplitude 



and the CP-asymmetry 



A, 



CP 



i \A\A+\^ + -m% \A-\^] 



r(E,^77)-r(p, ^77) 

T{Bs ^ 77) + r(5, ^ 77) 



(8) 



(9) 



we get 



D 



Acp = 2Im{\e 



(10) 
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where 





K) 


= aPi 




h) 


= aP2 + b 




-) 


= cP, 




-) 


= cP2 + d 



(11) 



and 



fBs A TT 9- 



\/27r " SniBAsiniB, - A 



{{mbA{mb) + msA{ms))D{fi) + m,A{mc)E{fi)} ] (12) 



In eq. (p^, we used the parametrization, 

C',ff{li) = Pi{li)\e + P2{li) (13) 

with 

A. = ^lel«el..|e^' (14) 
mtrrib 

Here, ^^^^j is chosen to be as real and i^^^b complex, namely ^^^^^ = \iN,bbW^ ■ Finally the 
functions Pi(/i) and P2{lj)^ in the LLog approximation |21|, are 



PM = v''^''F2{y) + ^{v''/''-v''^'')G2{y) 

mr 



3 

+g,(C5^°(/i) + N^c^'^iii)) + +Q^(!!^c^o^^^ ^ N^Uhc^o^^^^ 

+C2{mw)jZKr]''^ (15) 
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and 



D = |Ae|2{8|Ti(+^p + m|jjTi"^|2} + 2i?e(Ae){8i?e(TWT(+)*) 

+m%Re{Tt''T^'>)} + {SlT^+^l^ + m%\Tt^\'} (16) 



In eq. |15|, rj = a;s(mvi/)/as(yu) and hi, ai are numbers which appear during the evaluation of 
the Wilson coefficients p2| . 

Now, we would like to add the LD distance contributions due to the process -B^ ^ 07 — >■ 77 



TT[ . These effects can be taken into account by the redefinition of the functions Ti , T2 , T| , 



( + ) rp{ + ) rp{-) rp{-) 



rp'i + ) _ rpi + ) p ( + ) 

rp'i-) _ rp{ + ) p (-) 



Tt^ = T^ + AaS (17) 



where 



. _^^f.(0)/,(0)A, '"'''"^--'"^' 



,(+) 



^^-^ - v^^^F,(0)/,(0)A 



Note that there is also a contribution due to the chain process Bs —>■ (pip ^ </)7 ^ 77 and 
it is negligible compared to the one due to the decay -B^ — > 07 — > 77 
Finally, we derive the branching ratio for the given process as 

1 



(+) |2 + ^4 ^ I |2| + 4 „A rfi^(+) J,..(^i+h 



= — [|A,|^{8|rrf + m^^jTrf} + 4i?e(A,){8rri?e(T2 

64 7rmDrtot 



+m|, Tr^i?e(Tro} + {8|Trr + ^i^rr}] (19) 

In the 3ifDM(02)(see Appendix C and [|18|), the physical quantities under consideration 
are the same with the redefinition of Ag, 

A, = ^e^,,e^,,,(cos2 6 + t sin^ 9) (20) 
mtmb 

3 Discussion 

In the model III, there are many parameters, such as complex Yukawa couplings, ^[^j^ (i,j are 
flavor indices), masses of charged and neutral Higgs bosons and they should be restricted using 
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experimental results. All the Yukawa couplings except ^^^^ and ^^^^ are cancelled based on 
the experimental measurements by CLEO p^ , 



Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10"^ , (21) 



AF = 2 mixing and the p parameter |2^ (see pO[ for details). This discussion also allows us to 



neglect the contributions coming from the primed operator set. Further, the CP parameter "6" , 
appearing in the combination ^^^jC^bb = \^N,tt^N,bb\^~^^ y restricted due to the experimental 
upper limit on neutron electric dipole moment dn < 10~^^ e.cm, which gives an upper bound 
to the combination ^^Jrn{^N,tt^N*bb) < ~ ^00 Gev l2|. 

Now, we are ready to start with the discussion oi CP asymmetry in our process. Note that, 
in the analysis, jCy-^-^l is allowed to lie in the region, 0.257 < ICy^^l < 0.439 due to the CLEO 
measurement (see |^) and the parameters ^^^j, C^feb ^ restricted. We choose the scale 



as /i = ^ since we predict that this choice reproduce effectively the Next to Leading Order 



(NLO) result (see [TT 



In Fig |T](|), we plot the Aqp of the decay Bg 77 with respect to sin 6*, for ^^1,1, = 40 mf, 
and m//± = 400 GeV, in the case where the ratio \rtb\ = < 1, without (with) LD effects. 



■'N,bb 



in model III. Aqp is restricted in the region bounded by solid (dashed) lines for C^-^-^ > 
{Cj" < 0). Fi gures show that Acp is at the order of 10 ^, which is a weak effect. Acp is larger 
for Cj-^'^ > 0, as compared to the C^^^ < case and it is possible that Acp vanishes and even 
take negative values for sin 6' > and C^^^ < 0. Addition of LD effects enhances Acp small in 
amount. 

Fig ||(§) denotes the same dependence of Acp for 3HDM{02) with (without) LD effects. 
The area of the restricted region and the possible values of Aqp are smaller as compared to the 
2IIDM case. For example, for sin 6' = 0.5 and Cj^^ > 0, the upper limit of Acp decreases by an 
amount of 50%. However, the order of Acp still remains the same, namely 10~^. Furthermore, 
Acp is not sensitive to the charged Higgs boson mass (see Fig ^ and (^. 

Figures |7|-PD| show the Br of the given process for the model III and 3HDM{02). In Fig |^ 
(§), we present the dependence of Br on sin^ without(with) LD effects for \rtb\ < 1, in model 
III. The restricted region lies between solid (dashed) lines for C^^^ > {C^^^ < 0). Note that 
Br for C^-^-^ < is greater than the one for Cj^^ > 0. It is observed that the enhancement of Br 
in the SM is negligible (BrsM = 3.45 lO"'^ with LD effects and BrsM = 4.71 10"^ without LD 
effects). Therefore model III can not be distinguished from the SM with the measurement of Br 
of the given process, for \rtb\ < 1. Fig || and ([1^ denote the same dependence for 3HDM{02) 
and the results are similar to the 2HDM case. 
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Finally, in Fig 11, we plot the CP ratio i? = |^| with respect to sin 6' for 2HDM. Here, the 
solid line corresponds to Cj^^ > and the dashed line to Cj^^ < 0. This figure shows that R is 
not sensitive to CP violating parameter 6 and there is a small enhancement compared to the 
SM value {Rsm = 0.845). This ratio is also non-sensitive to the models under consideration. 

For completeness, we would like to note that there are some uncertainities coming from 
the choice of the bound state parameters nfj^ , and the decay constant fs^- The physical 
quantities are sensitive to these parameters. For example, the larger nfj^ (smaller A^), the 
larger Br, R and the smaller Acp- 

In conclusion, we study the CP asymmetry of 77 decay in the framework of the 

model III and 3HDM{02). Further, we analyze the Br and R ratio of the given process. We 
can summarize the main points of our results: 

• In model III, a weak Acp is possible and it is at the order of 10~^. This effect increases 
with the addition of LD contributions and this holds also in the 3HDM{02) model. The 
measurement of such a small value of Acp can give information about the sign of Cj^^ . 

• The Br ratio is not sensitive to CP violation parameter 6 and the enhancement as com- 
pared to SM is negligible in both models, for \rtb\ < 1 

• The R ratio is also non-sensitive to the parameter 6 in both models. 
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A Appendix 
Model III 



The Yukawa interaction for the 2HDM in the general case is 

where L and R denote chiral projections L{R) = 1/2(1 =p 75), (pi for i = 1,2, are the two scalar 
doublets, rjfj^ and ^fj^ are the matrices of the Yukawa couplings. With the choice of 0i and 



25 



1 

V2 



\ ,( V2x^ 



(23) 



and the vacuum expectation values, 



< 01 >= ^ ( ° j ; < 02 >= , (24) 
we can write the Flavor Changing (FC) part of the interaction as 

(^Y,FC = e{jQ^L02f/,/? + ^?,Q^L<P2D,R + h.C. , (25) 

where the couplings ^^'^ for the FC charged interactions are 

— ^neutral VcKM , 
^ch = VcKM ^neutral , (26) 

and C,neutrai Q defined by the expression 

= ivH'^r'^'^Vj^'^ . (27) 

Here, the charged couplings appear as a linear combinations of neutral couplings multiplied by 
VcKM matrix elements. 

B Appendix 
3HDM(02) 

In the 3HDM the general Yukawa interaction is, 

Cy = vl'jQiL'PiUjR + vSQiL(piDjR + ^I'jQiLhUjR + ^^.QnhDjR 

+ PijQiLhUjR + pf'jQiL<l)3DjR + h.c. , (28) 



^In all next discussion we denote £,neutrai ^n'^ ■ 
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where 0j for i = 1,2, 3, are three scalar doublets and vfj^y ^ij'^ ^ Pij^ Yukawa matrices 

having complex entries, in general. Now, we choose scalar Higgs doublets such that the first 
one describes only the SM part and last two carry the information about new physics beyond 
the SM: 



.0 



1 

71 



V2H+ 



1 

71 



V2F+ 



(29) 





V 



; < 02 >= ; < 03 >= 



with the vacuum expectation values. 

Note that, the similar choice was done in the literature for the general 2HDM (model III) 
The Yukawa interaction responsible for the Flavor Changing (FC) interactions is 

and the couplings ^^'^ and p^'^ for the charged FC interactions are 



(30) 



(31) 



p1 



VcKM , 
VcKM , 
PN VcKM 5 
yCKM PN , 



(32) 



and 



SAT 
U,D 

Pn 



L 7 S ^R 

U,D^-1 JJ^D y 



U,D 
R 



(33) 



In the 3IIDM model, the Higgs sector is extended and therefore the number of free pa- 
rameters, namely, masses of new charged and neutral Higgs particles, new Yukawa couplings, 
increases. Fortunately, by introducing a new transformation in the Higgs sector and taking the 
3HDM Lagrangian invariant under it, the number of free parameters can be reduced enormously 
T8| . Taking the following 0(2) transformation: 



cos a 02 + sin a 03 , 
—sin a 02 + cos a 03 



(34) 
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where a is the global parameter, which represents a rotation of the vectors 02 and 03 along the 
axis where 0i lies and assuming the invariance of the gauge and CP invariant Higgs potential 

^(01,02,03) = Ci(0+0i-t;V2)2 + C2(0+02)' 

+ C3(03+03)' + C4[(0|01 - t;V2) + 0^02 + 0^03]' 

+ C5[(0+0l)(02+02)-(0^02)(0^0l)] 

+ C6[(0+0l)(0+03)-(0^03)(03+0l)] 

+ C7[(0+02)(0^03)-(0^03)(0^02)] 

+ C8[i?e(0+02)]' + C9[/?e(0+03)]' + cio[i?e(0+03)]' 

+ cii[/m(0+02)]' + ci2[/m(0+03)]2 + ci3[/m(0+03)]2 + ci4 (35) 

we get the masses of new particles as 

mp± = mii± , 
mjy4 = mj^2 , 

(36) 

Further, the application of this transformation to the Yukawa Lagrangian eg . (^81) allows us to 
write the equality 

where 

clU{D) (.U{D) . U{D) ■ 

Kij = ^ij ' COS a + 'sma , 

tU{D) ^U{D) . . U{D) /ooN 

Pi/ = 'sma + Pi/ 'cos a. (38) 

and therefore the Yukawa matrices and p^(^) can be parametrized as, 

= e^sinO , 

p^ = ie^sine. (39) 
Here e^^^^ are real matrices satisfy the equation 
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and T denotes transpose operation. Finally, we could reduce the number of the Yukawa matrices 
^^'^^^ and p^'^^\ by connecting them with the expression given in eg . (p9D . Further, we take 
into account only the Yukawa couplings ^^^t, ^Nbb^ P^tt PNbbi since we assume that the 



others are small due to the discussion given in |20 



C Appendix 

The operator basis and the Wilson coefficients for the 
decay h S77 in the model III 



The operator basis is the same as the one used for the 6 — > 37 decay in the model III pO| and 
SU{2)l X SU{2)r X f/(l) extensions of the SM 



01 = {sLal^,CL^3){cL^3l^hLa), 

O3 = (sLalt^bLa) J2 i^LflYqLfs), 
q=u,d,s,c,b 

q=u,d,s,c,b 

q=u,d,s,c,b 

q=u,d,s,c,b 

g 

^7 = -—^Sa(J^^{mbR + msL)baT^"', 

Oio = {sLali,CLa){cR(3l^hRp), (41) 

and the second operator set 0[ — 0[q which are flipped chirality partners of Oi — Oiq: 

O'l = isRc,'yf,CRf3){cRf3YbRa), 

02 = isRa-ff^CRa)icRf3YbRI3), 

03 = (sRaif^bRa) {qRPi^qRP)^ 

q=u,d,s,c,b 

04 = {sRa^^hRp) Y i^Rpl'^QRa), 

q=u,d,s,c,b 

q=u,d,s,c,b 
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q=u,d,s,c,b 

g 

O'g = {sRa'l^CRp){cLpYbLa) , 
Oio = {sRa-i^CR^){cLpYhLI3) , (42) 

where a and P are ^[/(S) colour indices and J-'^'^ and ^^'^ are the field strength tensors of 
the electromagnetic and strong interactions, respectively. In the calculations, we take only the 
charged Higgs contributions into account and neglect the effects of neutral Higgs bosons (see 
||T6| for details). Further, in our expressions we use the redefinition. 



e^'^ = ^^e^'^. (43) 

Denoting the Wilson coefficients for the SM with C^'^{mw) and the additional charged 
Higgs contribution with C^{mw)^ we have the initial values for the first set of operators 
(eq.( ^)) (1^ and references within) 

<^1^3!..6,9,lo('^Vl/) = , 

Cl'\mw) = 1 , 



4(x - 1)4 24(x - i; 



3 



r^SM I \ "-"-^ 1 X -t- OX -t- 

^« = ~4(^^^""+ 8(x-l)3 ' 

CfJ. . 6,9,10 ("^W) = , 

^ts ^tb 
^t^b yts ^tb 

^ts ^tb 

+ -V + ^^"^-^^ + ^N.b'^^G.iy) , (44) 

^t^b yts ^tb 



and for the second set of operators eq. ( ^2]) 
<^Wo(^iy) = 0, 

^?..6,9,10(^W') = 0, 

C.^'imw) = —A^lbs^.+^lss){^lbb + ilsb^)Fi{y) 

vts vtb 
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mtrrib Vu ' Vth 



mf ' Vts ' Vtb 



+ ^ (^"^.^^^ + + , (45) 

TTT't^b Vts ^tb 

where x = mf/m^r and y = mf/mjj±. The functions Fi{y), F2{y), Gi{y) and 6*2(1/) are given 
as 

72(1/ -1)3 12(y-l)4 ^' 

^(%-3) y(-% + 2) 
^'^y^ = 12(^^+ 6(y-l)3 

^^^^ 24(?/-l)3 4(y-l)4 ^' 

Note that we neglect the contributions of the internal u and c quarks compared to one due to 
the internal t quark. In the 3HDM{02) model, the Wilson coefficients are obtained with the 
replacement ^^^^^ — > e^(-^) and the redefinition of Xg, \g = m/m^ ^N,tt^N,bbi'^'^^'^ ^ + isin^O). 
For the initial values of the Wilson coefficients in the model III (eqs. (|^)and (^)), we have 

^i,f^.6!9,io("^M^) = ' 
Cr^'\m^) = 1, 

<-'l,2,3,...6,9,10l"^W^J - U , 



Cf^^^(mH.) = Cf*^(mH.) + Cf (mvK) , 

C'i''''^'{mw) = C'i'\mw) + C'i'{mw) . (47) 



At this stage it is possible to obtain the result for model II, in the approximation ^ 
and ^ ~ 0, by making the following replacements in the Wilson coefficients: 

o 1 



^st*^Fb = -mtm,, (48) 



and taking zero for the coefficients of the flipped operator set, i.e C- — > 0. 
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The evaluation of the Wilson coefficients are done by using the initial values (^2-ffDM(3/i'Z)Af(02) 
(Cj'^^'^^^''^^'^*^'''^^^'') and their contributions at any lower scale can be calculated as in the SM 
case 
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D Appendix 

The necessary functions used in the Wilson coeffi- 
cients 

The explicit forms of the functions /(m^), J{mq) and A(mq) appearing in eqs. ^ and ^ are 



2 

TfL 

liniq) = l + ^A(m, 



m 



Bs 



- Ami 





Ami 

) 


\/^k - 


Ami 


J Ami- 


-^k 



^k 



me, - Jm%^ - Amj, / ^m^ 



A(m„) = - I 2arctan(-^^ —) - tt 1 for ^ < 1. (49) 

ms, I Ami 
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Figure 1: Acp as a function of sin 9 for ^^55 = 40 and mjj = 400 GeV in the region \rti,\ < 1, 
at the scale fi = 171^/2, without LD effects, in model III. Here Acp is restricted in the region 
bounded by solid lines for C^-^'^ > and by dashed lines for C^^-^ < 0. 



18 




19 




0.1 0.2 0.3 0.4 0.5 

si7i 9 



Figure 4: The same as Fig. ^ but with LD effects. 
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Figure 5: Acp as a function of mH± for sin 6 = 0.5 and ^^^^ = 40 in the region \rtb\ < 1, at 
the scale /i = mb/2, with LD effects, in model III. Here Acp is restricted in the region bounded 
by solid lines for C^^^ > and by dashed lines for C^^^ < 0. 
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Figure 6: The same as Fig ??, but for 3HDM{02) and 
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Figure 7: Br as a function of sin 9 for = 40 rrii, and = 400 GeV in the region |rt6| < 1, 
at the scale = without LD effects, in model 111. Here Br is restricted in the region 

bounded by solid hnes for Cj-^^ > and by dashed lines for C^^-^ < 0. 
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Figure 8: The same as Fig |^ but with LD effects. 



4.72 



4.715 - 

4.71 
4.705 - 
4.7 - 
4.695 - 

4.69 
4.685 ^ 

4.68 - 
4.675 - 

4.67 







0.1 



0.2 0.3 

sin 6 



0.4 



0.5 



Figure 9: The same as Fig but for 3HDM{02) and e^ ;,^ = 40 m^. 
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Figure 10: The same as Fig but with LD effects. 
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Figure 11: R ratio as a function of sinO for ^^^^ = 40 m;, and = 400 GeV^ in the region 
\rtb\ < 1, at the scale /i = mf,/2, in model III. Here R ratio is restricted on the solid (dashed) 
line for C^^^ > {Cj^^ < 0). Note that, for 3HDM{02), R ratio is almost the same as the one 
calculated in model III. 
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